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ABSTRACT 



An Observing System Simulation Experiment (OSSE) is 
conducted to investigate the sensitivity of the Naval Research 
Laboratory (NRL) mesoscale model to atmospheric water and 
surface wind speed input. The explosively deepening ERICA 
IOP-4 storm which occurred off the North Carolina coast on 4 
January 1989 is used as a case study. Wind speed and water 
vapor output from the full physics version of the model is 
shown to reasonably simulate data available from the Special 
Sensor Mi crowave/ Imager . Assimilating this modeled satellite 
information into the degraded version of the model shows ERICA 
IOP-4 to be a dynamically driven storm. 

Assimilation of simulated wind speed observations yields 
the best performance for the degraded model in terms of 
surface pressure deepening. The model exhibits sensitivities 
to the coverage area of updated information. Two treatments 
for atmospheric water are tested, showing the sensitivity of 
the model to these inputs despite their secondary importance 
in forecasts of storm deepening. 
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INTRODUCTION 



A. MOTIVATION 

To improve numerical weather prediction, satellite data 
needs to be assimilated with other model inputs. Efficiently 
incorporating these remote observations into numerical weather 
prediction models improves forecasts by bridging the gap 
between simulation and the real world. Problems arise in 
marrying a sensor of fixed temporal and spatial resolution to 
a model whose stability requirements dictate a specific finite 
resolution scheme. The objective of this study is to examine 
the potential for incorporation of DMSP Special Sensor 
Microwave/ Imager (SSM/I) data into the Naval Research 
Laboratory's Limited Area Dynamical Weather Prediction Model. 
Products to be tested include SSM/I derived wind speed and 
atmospheric water content. The case study to be examined is 
an explosive storm occurring on 4 January 19 89, part of the 
Experiment on Rapidly Intensifying Cyclones over the Atlantic, 
Intensive Observation Period Four (ERICA IOP-4) . Operational 
models failed to predict this storm's extreme deepening rate 
of 60 mb in 24 hours. By nudging the modeled development of 
this storm with information derived from the simulation of 
three SSM/I passes, a prediction closer to the final intensity 
of the storm, valid at 0000 UTC 5 January 1989, may be obtained. 
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Post-storm analysis (Sanders, 1990) places the 0000 UTC 5 
January pressure of the storm at 936 mb, while the National 
Meteorological Center (NMC) analysis gives a value of 956 mb. 
A disparity evolves between the predictions and the actual 
storm as development progresses. The widening gap motivates 
inquiry into the most important factors in the storm' s rapid 
development, and the possibility of effectively simulating 
satellite observations. 

The SSM/I is an excellent candidate for providing observed 
values for validation and improvement of model predictions. 
Fully operational, the SSM/I provides reliable atmospheric 
water parameters suitable for operational forecasting 
(Alishouse et al . , 1990a) . The improved non-linear algorithm 
has been validated by radiosonde observations, and shows 
marked sensitivity to cloud liquid water content, rain and sea 
ice. Cloud liquid water has been validated by the NOAA-WPL 
profiler network (Alishouse et al . , 1990b). While not 
reliable over land, the derived algorithm is highly effective 
at determining marine cloud liquid water. For the North 
Atlantic Ocean domain used in this study, the data provided is 
suitable. The SSM/I also provides surface wind speed over the 
ocean based on wave structure and foam coverage. The 
information is usable in the modeling case, but wind speed 
data is not as accurate as other parameters (Hollinger, 1991) . 
Sources of SSM/I error will be discussed in Chapter II. 
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The SSM/I's ability to provide global, day and night 
observations make it a potentially important model input. Use 
of microwave imagery is advantageous since microwaves 
penetrate cloud coverage. Microwaves are independent of solar 
illumination, and are readily comparable with visible and 
infrared images (Hollinger and Sandlin, 1990) . Microwave 
transmission is sensitive to the physical state of water in 
the atmosphere, dependent upon humidity, clouds and rain. 
Proven operational use of SSM/I data includes the marine 
weather updates covering the approach of Hurricane Hugo, which 
incorporated 3 -hour old SSM/I wind speeds. Fleet Numerical 
Oceanography Center (FNOC) also produces SSM/I derived wind 
speed fields. The Naval Research Laboratory -West (NRL-West, 
formerly the Naval Environmental Prediction and Research 
Facility, NEPRF) has concluded that SSM/I wind speeds are 
equivalent in quality to ship reported winds, and is now 
attempting to successfully assimilate these winds into NOGAPS 
(Goerss, 1988) . Data collection has begun with the SSM/I 
follow on platform, the SSM/T2 atmospheric sounder. 

B. BACKGROUND -EXISTING THEORIES OF EXPLOSIVE CYCLOGENESIS 

While not a comprehensive review of the subject, several 
theories on aspects of cyclogenesis that can provide the 
impetus for an explosive event are mentioned here to provide 
direction in the post experiment analysis. Fosdick and Smith 
(1991) have listed processes previously identified as 
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important to explosive cyclones, or "bombs.” These include 
upper level cyclonic vorticity advection, upper level 
divergence associated with a well defined jet streak, boundary 
layer fluxes of heat and moisture, latent heat release and low 
static stability as a precursor to development. Their 
synoptic discussion of an explosive storm developing in 
roughly the same area is considered here. The initial setting 
consisted of a surface low in a strong baroclinic zone 
downstream of a large amplitude upper level trough, near the 
entrance region of a jet maximum. Midway through the 
explosive phase, strong low-level thermal advection was 
occurring within the cyclonic circulation. By the end of the 
explosive phase, this thermal advection had extended southwest 
and northeast of the center. 

The precipitation pattern of the Fosdick and Smith storm 
changed according to the phase of the storm. At the beginning 
of the explosive phase, precipitation was predominantly 
convective and located in the warm sector. At the midpoint of 
the explosive phase, amount and areal coverage of 
precipitation had increased. Precipitation was now a maximum 
over the low and still mostly convective. Relating the amount 
of rainfall to latent heat release, Fosdick and Smith were 
able to draw conclusions about the importance of latent heat 
release on explosive cyclogenesis. Latent heat release 
contributed to height falls below the level of maximum 
heating. Latent heat release accounted for 33% of total 
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development at 700 mb. Lowering the level of maximum heating 
increases cyclone development. 

Studies of interaction between baroclinic and diabatic 
processes suggest that a nonlinear relationship exists between 
latent heat release and rapid development (Kuo et. al, 1991) . 
Kuo et al . ' s study of the QEII storm showed the surface 
cyclone related to quasi -geostrophic forcing of a mid- 
tropospheric short wave trough for periods of greater than 12 
hours. Development was also associated with intense mid-level 
and upper- level vorticity maxima. Moist f rontogenesis 
contributed to rapid marine cyclogenesis. The conclusion of 
Kuo et al . ' s study was that rapid cyclogenesis must be viewed 
in the context of moist baroclinic instability with strong 
nonlinear interaction between dynamic and diabatic processes. 

Varying the treatment of atmospheric moisture can have a 
distinct impact on model predictions. In a study of nine 
Western Atlantic cyclones, the Arakawa- Schubert scheme gave 
deeper, more realistic deepening rates than the Kuo convective 
parameterization scheme (Kuo and Low-Nam, 1990) . 
Precipitation is driven by mesoscale slantwise ascent in the 
vicinity of the warm front in a rapidly deepening cyclone. In 
concurrence with the previously discussed studies, latent heat 
release (LHR, here associated with the sloping frontal cloud 
band) promotes f rontogenesis and increased low-level vorticity 
when LHR interacts with adiabatic frontal circulation. Models 
which can determine LHR from internal dynamics should give 
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better results than models using arbitrarily prescribed 
parameterization schemes (Kuo and Low-Nam, 1990) . Thus, using 
simulated observations of precipitation to drive latent 
heating release rates should improve model simulations. A 
further result of the Kuo study was that significant amounts 
of water vapor were needed in the initial stages to support 
rapid deepening. 

Kuo and Low-Nam's conclusion that dynamically formed 
cyclones are more easily predicted than diabatically formed 
storms is supported by studies based on NMC's Nested Grid 
Model. Sensitivity tests using the NGM Model examining the 
role of latent heat release in IOP-4 attributed 27% of its 
deepening to latent heat release, suggesting that the storm 
was primarily driven by dynamics probably associated with the 
strong short wave aloft (Pauley et. al, 1991). 

C. SYNOPTIC DISCUSSION OF THE ERICA IOP-4 CASE STUDY 

The ERICA IOP-4 storm of 4-5 January 1989 was the deepest 
cyclone occurring during the ERICA project (Hartnett et . al, 
19 89) . Between 09 00 and 1500 UTC (all times given are 4 
January 1989 unless specified) the storm deepened at the 
explosive rate of 24 mb in 6 hours, enroute to a final 
estimated pressure of 936 mb. Figure 1 presents the storm's 
deepening rate and low pressure center track. The IOP-4 storm 
was the trailer of a pair of storms depicted in Figure 2. 
Rapid 
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deepening occurred immediately upon upper air forcing reaching 
the coastline at 0000 UTC. Strong mesoscale development in 
the cyclone interior during the explosive phase resulted in a 
tight cloud band surrounding the clear center. Satellite 
imagery implied that development began in the northwest 
quadrant of the cyclone around several equally strong low 
pressure centers at 1200 UTC. These centers made a partial 
revolution around the previous cyclone center. The 
intensification of this low coincided with the period of most 
rapid deepening. 

A detailed discussion of synoptic and mesoscale 
characteristics of the storm lays the groundwork for 
comparison of model experiments. Evolution of sea level 
pressure deepening began with a 996 mb low east of Cape 
Hatteras, North Carolina (Nieman et al . , 1991). Wann and 
moist air over the Gulf Stream was separated from cooler and 
drier continental air over the coastal plain by a coastal 
trough extending southwestward from the cyclone center (Figure 
2) . Intensification began as a surface cold front, coupled 
with an upper level jet -front system, advanced from the west. 
Strong southwesterly winds developed east of the coastal 
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Figure 2. NMC analyzed surface pressure (solid lines, mb) 
and 1000-500 mb thickness (dashed, dm) at 0000 UTC 4 January 
1989. 
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trough in the Gulf Stream boundary layer due to mixing of 
higher momentum aloft downward. These strong southwesterly 
winds combined with negative air- sea temperature differences 
contributed to combined sensible and latent heat fluxes of 
1500 W/m^ upwards from the ocean surface. Ahead of the cold 
front, a deep convective cloud line developed over the upward 
heat fluxes. The resulting flow transported warm moist air 
northward toward the warm f rontogenesis region. Cyclogenesis 
commenced beneath the 500 mb flow and strong 700 mb ascent 
{Nieman et. al, 1991). 

Due to the linear pressure drop shown in Figure 1, Niemann 
et. al could not identify a prominent deepening phase. In the 
first 6 hours of the case study, the intensifying cyclone was 
centered near the triple point, while an incipient bent-back 
frontal feature was developing to the west. The mesoscale 
cyclone continued to trail farther and farther behind the 
triple point and deeper into the polar air stream. The 
cyclone exhibited a scale and structure similar to polar low 
type cyclogenesis as it propagated from the southwest in a 
reverse shear environment. 

Following temperature evolution at 850 mb, Neiman et. al . 
were able to identify four phases of frontal evolution 
developed by Shapiro and Keyser (1990) . Specifically; 

1. Broad Baroclinic phase 

2 . Frontal Fracture phase 
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3. T-bone phase, and 

4. Warm Seclusion phase 

The broad east -west baroclinic zone present at 0000 UTC 
supported incipient cyclogenesis. By 0600 UTC progress of the 
second phase was evident as the warm and cold fronts were 
distinctly fractured at the triple point. At 1200 UTC, a T- 
bone frontal configuration developed consisting of the north- 
south oriented cold front intersecting the east -west warm 
front. Finally, between 1800 UTC and 0000 UTC 5 January, the 
cyclone attained its lowest pressure and maximum vorticity in 
the warm core seclusion phase. Cold polar air nearly 
encircled the cyclone center by 0000 UTC (5 January) , 
secluding a warm pocket near the low center. At the end of 
the period, the outward sloping bent -back occlusion was 
distinguishable through 500 mb, although it weakened with 
height (Niemann et. al, 1991) . 

Flight level observations documented mesoscale features, 
producing temperature and streamline data. Based on cross- 
sectional analysis, the 0600 UTC cold front was nearly 
vertical below 800 mb. Cold front convective bands contained 
20 km scale eddies and periodic gaps in the along- front 
direction. The warm front consisted of a steep leading edge 
with a positive isentropic potential vorticity anomaly at the 
leading edge of the front. Radar cross-sections revealed a 
shallow slope along the warm front while under stratiform 
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precipitation; and a nearly vertical slope in the vicinity of 
the mesoscale convective elements. Vertical portions of the 
front displayed enhanced cross- frontal convergence. 

Streamline analysis from the same aircraft data shows 
strong confluence along the warm convective band, with weaker 
confluence aligned with the cold frontal band. The primary 
cyclone center at 0600 UTC was located in the closed 
circulation approximately 40 km west of the triple point. A 
secondary circulation center existed 125 km west of the triple 
point along the bent -back front. The analysis clearly 
displays multiple vortices during the early stage of 
cyclogenesis. One of the objectives of this study is the 
investigation of how assimilation of satellite data can 
improve the model's description of these features. 
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II. TREATMENT OF SATELLITE DATA 



A. SENSOR CONFIGURATION AND ALGORITHMS 

The Special Sensor Microwave /Imager flies on board the 
Defense Meteorological Satellite Program (DMSP) block 5D-2 
spacecraft F8, operational since 1987. The spacecraft orbits 
at 833 km in a sun- synchronous , near polar orbit resulting in 
a 1400 km swath width. The sensor consists of a seven 
channel, four frequency (19.35, 22.235, 37.0 and 85.5 GHz), 
linearly polarized passive microwave radiometric system. 
Continuous rotation at 31.6 RPM about a vertical axis allows 
the sensor to measure brightness temperatures over a 102.4° 
angular sector. The imaging scanner operates left to right. 
During each scan, 128 radiometric samples are taken at the two 
85.5 GHz channels, with 64 samples taken on alternating scans 
at the remaining five lower frequency channels (Hollinger et 
al., 1987). 

Scan geometry will influence the nature of geolocation 
error in an SSM/I image. The satellite track moves along the 
Y-direction at 6.58 km/sec, rotating about the vertical as 
described. This results in a separation distance between 
successive scans of 12.5 km along the track, which is nearly 
the resolution of the 85.5 GHz channels. Using a sampling 
interval of 4.22 msec (the time for the beam to travel 12.5 km 
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across the track) increases the along -scan beam diameter. As 
a result, the 85.5 GHz beam is nearly circular. The lower 
frequency channels sample at a 8.44 msec interval on 
alternating scans. Figure 3 depicts the footprint of beams in 
a swath, indicating that beam overlap occurs predominantly at 
the edge of the swath (Hollinger et al . , 1987). 

Beam inefficiency coupled with scan geometry contributes 
to SSM/I error. The most efficient beams (beams having the 
highest concentration of energy in a 3 dB bandwidth) are the 
19.35 GHz channels. The vertical beam is 96.1% efficient, and 
the horizontal is 96.5% efficient. Least efficient is the 
85.5 GHz horizontal, calculated at 91.1% based on the 
deployment of antenna components, alignment of the sensor with 
the spacecraft, and feedhorn design. NRL computations have 
determined an along beam pointing accuracy of 5.3 km, and a 
cross track accuracy of 3.2 km, both less than the 85.5 GHz 
beam half width (Hollinger et al . , 1987). 

Noise Equivalent Temperature Differential (NEaT) 
represents the radiometer sensitivity. To reduce radiometer 
gain fluctuation and, subsequently, NEaT, the instrument 
undergoes a radiometric calibration every 1.19 seconds. 
Calibration is done by passing the feedhorn beneath two fixed 
calibration reference targets, a 300K hot receiver and a 3K 
cold cosmic background view. Scene incident temperature is a 
linear function of these calibration readings and received 
brightness temperatures. Five samples from each scan are 



14 



Along Track 
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3-db Beam Diameter 
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Figure 3. SSM/I 37 Ghz and 85 Ghz channels beam 
footprints within a swath for various angles from the 
nadir (Hollinger et al . , 1987). 
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averaged to reduce noise estimates in hot and cold load 
readings. Calibration parameters are updated on each scan to 
compute scene brightness temperatures. The accuracy of scene 
brightness temperatures incident on the reflector can improve 
with several corrections applied. These corrections account 
for feedhorn spillover, cross polarization coupling and 
sidelobe contribution. 

1. Wind Speed Algorithm 

The recommended surface wind speed algorithm is 
linear, with varying degrees of accuracy (Hollinger, 1991) : 



SW=1^1 .90+1.0969* .455* ^ 6 * .7860* 

where SW is wind speed in meters per second. Typical wind 
speed errors are near 2 m/s. Tbj 9 v# are the 

brightness temperatures in the 19 GHz, 22 GHz and 37 GHz 
vertically polarized beams. T 37 H is the brightness temperature 
in the 37 GHz horizontally polarized beam. 

Microwave radiation, particularly at the lowest three 
frequencies of the SSM/I, is heavily attenuated by rain, which 
in turn dampens the wave -generated wind speed signature. 
Table 1 lists rain rate dependent wind speed accuracies . 
Accuracy of the wind speed retrieval drops off severely with 
rain due to sensor frequency limitations, not algorithm 
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TABLE 1. WIND SPEED ACCURACIES BASED ON RAIN FLAGS 
(Hollinger et al . , 1989) 



Rain Flag 


Criteria 




Accuracy 


0 


T 37 VT 37 H > 
AND 
Ti9H < 


50 


< 2 m/s 


1 


T 37 V-T 37 H < 

OR 

Ti9h > 165 


50 


2-5 m/s 


2 


T 37 V-T 37 H < 


37 


5-10 m/s 


3 


T 37 V-T 37 H < 


30 


> 10 m/s 



problems. (Hollinger, 1989) 

Figures 4, 5, and 6 show the rain flag distribution and 
geographic position of the three swaths. The most severe wind 
speed error in this case study is 5 m/s. 

2. Water Vapor Algorithm 

The recommended water vapor (total precipitable water) 
retrieval algorithm is non-linear in 22 GHz (Hollinger, 1991) : 

m^=232 . 89393 -.14 8596 829125 *T 22 /. 006193* (^ 22 ^) 2 

-3 .6954*r„ 

where WV is the total precipitable water in kilograms per 
square meter. This algorithm results from regression of 
international radiosonde network data. Typical errors in 
total precipitable water are within 2.4 kg/m^. 
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Figure 4. Wind speed accuracies based on SSM/I rain flags 
for the 0000 UTC swath. 
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Figure 5. Wind speed accuracies as in Figure 4, but for 
the 0930 UTC swath. 
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Figure 6. Wind speed accuracies as in Figure 4, but for 
the 2200 UTC swath. 
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The validation report (Hollinger, 1991) recommends 
employment of a precipitation screen prior to computing water 
vapor over the ocean. Application of this screen obliterates 
most of the image, while neglecting this screen gives useful 
structural information. In the unscreened image, however, the 
magnitudes of precipitable water are much higher than the NRL 
model predictions, as discussed in Chapter IV. 

B. ERROR SOURCES IN THE SSM/I DATA 

The SSM/I was designed to meet an operational spatial 
resolution requirement of 14 km. However, the sensor has not 
always performed to this standard, displaying errors up to 30 
km near coastlines. The sensor's geolocation algorithm 
assumes that the spacecraft is always perpendicular to the 
earth- satellite vertical, and that the earth is a perfect 
oblate spheroid. Errors can be induced by many elements, 
including : 

• Spacecraft orbital elements 

• Spacecraft ephemeris data 

• Sensor pixel location algorithm 

• Sensor deployment /alignment of spin axis 

• Sensor alignment to spacecraft 

• Spacecraft altitude 
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Also included are errors in maps to which the satellite data 
must be compared (Poe and Conway, 1990) . 

Orbital element errors, errors in the set of parameters 
defining the size, shape and orientation of the orbit, are 
small contributors to SSM/I geolocation error. Comparisons of 
NORAD and NAVSPASUR data indicate the orbital elements 
contribution to be less than 1.5 km. (Poe and Conway, 1990) . 
Errors in the satellites tabulated position (ephemeris) are 
less than 3 km as computed by contrasting FNOC and NORAD 
models. Comparisons of FNOC and NORAD data displayed a 13 km 
subsatellite position difference coincident with FNOC^s 
dataset update time. This indicates the scale of the error 
induced by the data collection process. Further comparisons 
with ephemeris data from the U.S. Air Force 1000th Space 
Operations Group (1000th SOG) indicate differences between 
spacecraft and FNOC ephemeris can be as large as 15 km 
(Hollinger, 1991) . 

The geolocation algorithm employs many numerical 
approximations to reduce the computer processing load. This 
induces a small error in the satellite's tabulated position. 
Total error from this source does not exceed 4 km. 

The final NRL Calibration/Validation study of SSM/I data 
used the 85 GHz channel to study geolocation error in SSM/I 
derived images. This channel provides the best spatial 
resolution, 12 km. The largest geolocation error offsets the 
data 20 to 30 km in the eastern Mediterranean Sea, with common 
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smaller errors (15-20 km) existing near coastlines in the Red 
Sea. Computational methods have been proposed to reduce the 
errors to 1 to 3 km, but these processes have not been 
employed pending further testing with more extensive data 
sets. (Hollinger, 1991) . Table 2 summarizes SSM/I error 
estimates . 

TABLE 2. SSM/I PARAMETER ERRORS 



SSM/I ERROR 


CHANNEL RESOLUTION 


Orbital Elements 1 km 


85 GHz 12 km 


Ephemeris 15 km 


37 GHz 27.7 km 


Approximations 4 km 


22 GHz 46.5 km 


Total Geolocation 20 km 

error 


19 GHz 54 km 


SW/TPW algorithm error 


20 -H 27.7 = 47.7 km 
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III. MODEL DESCRIPTION 



A. DESCRIPTION OF NRL MODEL 

The NRL Limited Area Dynamical Weather Prediction Model 
accurately depicts processes and characteristics of the lower 
troposphere (Madala et al . , 1987). Using quasi -hydrostatic 
sigma coordinates permits reasonable modeling of baroclinic 
development. Other physical processes, notably cumulus 
convective heating, are parameterized. A brief description of 
the model provides a framework for understanding the structure 
of the sensitivity tests, which are described in Chapter IV. 

The model employs the Arakawa C-grid since that scheme 
presents several physical advantages (Arakawa and Lamb, 1977) . 
The C-grid best simulates the geostrophic adjustment process, 
which influences how the primitive equation system responds to 
ageostrophic initial conditions. When initial conditions 
specify an imbalanced region, gravity waves propagate out of 
that anomalous region. Use of the C-grid also conserves the 
integral properties of the atmosphere. Staggering the grid 
permits the separation of various waves by phase speed. 

At specified mass points on the grid, the quantities 
temperature, geopotential, specific humidity, and sigma are 
computed. U- velocity (east -west) is computed midway between 
latitudinal grid points, while v- velocity (north- south) is 
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computed midway between longitudinal grid points. The seven 
governing equations employed in the model are: 

• U- momentum 

• V- momentum 

• Thermodynamic equation 

• Moisture continuity equation 

• Surface Pressure Tendency 

• Hydrostatic equation 

• Continuity equation 

Integration of these seven equations over time by a split - 
explicit method generates the state variables. Use of the 
split -explicit scheme provides for separation of components 
into Rossby and Gravity wave modes. The primitive equations 
must be written in matrix form to achieve an eigenvector 
solution whose elements will satisfy the various propagation 
speed requirements (Madala et al . , 1987). 

B . MODEL STRUCTURE 

The main program calls several subroutines, each of which 
computes a state variable, a correction term, or paramaterizes 
a physical process. The objective of the program is the 
forecasting of meridional wind, zonal wind, temperature, 
surface pressure, geopotential height and sigma for the five 
prognostic and two diagnostic equations. From the 

initialization state chosen by the user, the model reads 
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fields of the wind components, air temperature, sea surface 
temperature, moisture, sea level pressure and boundary layer 
parameters. The eigenvector solutions that result from 
solving matrix forms of the primitive equations provide 
pressure gradient and divergence correction terms. Separate 
subroutines compute mass divergence, geopotential corrections, 
diffusion and Monin-Obukov boundary layer parameters. 

Of particular interest in this study, the subroutine 
CUPARA determines cumulus convection based on a modified Kuo 
scheme. According to this algorithm, convection exists only 
when the total moisture required for cloud formation is 
positive (Kuo, 1974) . This depends on cloud -environment 
temperature and moisture differences. New environmental 
values of temperature and humidity are fed back to the main 
program. In the absence of convective heating, stable heating 
is indicated if the relative humidity is greater than 96% 
(Holt et al, 1990) . Section D describes in greater detail the 
flow of the cumulus parameterization. 

C. TIME INTEGRATION 

Gradient processes in the atmosphere potentially trigger 
gravity wave disturbances. These disturbances are low 
amplitude waves with little energy that require a small time 
step if they are to be captured by the model. Deviation terms 
in divergence and pressure tendency contain elements that vary 
over the time scale of Rossby waves and require a larger time 
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step for proper modeling. The split -explicit method handles 
both wave scales by dividing the time step into sub- intervals 
and carrying out explicit integration within each sub- 
interval. Employment of this scheme in a three dimensional 
model in grid point space runs into practical difficulties 
(Madala et al . , 1987). 

On the grid, the array of gravity waves modes satisfy an 
array of CFL criteria. Each physically occurring gravity wave 
requires a different time step, as do the natural modes 
generated by the numerical model (one at each level) . 
Fortunately, the natural gravity modes form a complete set of 
eigensolutions to the boundary conditions. This set of 
solutions can be accounted for with constants, allowing 
variation of the dependent variables to be expressed as linear 
combinations of the structure functions (Madala et al . , 1987) . 

D. CUMULUS PARAMETERIZATION 

The cumulus parameterization subroutine defines convection 
and adjusts the state variables of temperature and moisture 
where cumulus convection exists. For this subroutine the 
model is strongly dependent on the initial magnitudes and 
distributions of the "observed" state values. In CUPARA, 
weighting functions for vertical heating as well as moisture 
divergence are computed for each model level. After vertical 
integration, the derived total moisture convergence 
(divergence) is added to any supersaturation that may have 
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occurred in the lowest level of the modeled atmosphere. 
Saturation vapor pressure, saturation specific humidity and 
relative humidity are derived for each level, as well as 
temperature and specific humidity at points of upward vertical 
velocity. Comparison of all levels yields the vertical 
gradient of moist stable energy. 

If the required amount of moisture to form clouds exists 
and the atmosphere is unstable, convective heating and clouds 
are indicated. When relative humidity exceeds 95%, no 
additional moisture is placed into the surrounding 
environment. Finally, the subroutine computes the change in 
temperature and humidity due to the mixing of clouds in the 
environment (Madala et al . , 1987). 

E. MODIFICATIONS TO THE MODEL 

For this study the model domain focuses on the western 
Atlantic Ocean, bounded by 20°-55° N and 40°-80° W as shown in 
Figure 7. The numerical grid consists of a 136 by 141 point 
horizontal plane (1/3 degree latitude by 1/4 degree 
longitude) , with 10 (version MIO) or 16 (version EE16) 
vertical levels. All experiments are initialized by NMC 
Hemispheric Analysis valid for 0000 UTC 04 JAN 1989. This 
series incorporates non- divergent winds on pressure surfaces 
and observed temperatures . Boundary values are nudged towards 
the 12 hour NMC analysis using: 
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=(l-a) 




+Y 

m 






where X is a dependent variable, subscript 'm' denotes the 
model tendency, r is a relaxation parameter, and represents 
the large scale (initialization field) value. Besides the 
cumulus parameterization scheme already discussed, 
parameterized processes include convective precipitation, dry 
convective adjustment, and Planetary Boundary layer (PBL) 
mixing (Holt et al . , 1990). 

Two versions of the model are used in this experiment 
series. The degraded MIO model provides the reference 
atmosphere. The bottom layer encloses the PBL, parameterized 
by the bulk aerodynamic (drag coefficient) method. Sea 
surface temperatures are from January climatology on a 1® 
resolution grid (Reynolds, 1982) . The sixteen layer EE16 
version of the model generates the simulated SSM/I 
observations. EE16 models the PBL using turbulent kinetic 
energy (TKE) closure theory, incorporating prognostic 
equations for turbulent kinetic energy transport and 
dissipation. TKE is derived as a function of shear 
production, buoyancy production and turbulent transport. The 
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surface energy balance is maintained via the soil slab model 
{Chang, 1979) . Dissipation of the energy results from 
imbalances in production, destruction and transport (Holt et 
al. , 1990) . 

Tests of MIO involve the assimilation of moisture and wind 
fields from EE16. Varying adjustments of assimilation time, 
spatial domain (simulating satellite swath widths) and the 
vertical distribution method for moisture comprise the 
sensitivity tests, described in the next chapter. 
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IV. OBSERVING SYSTEM SIMULATION EXPERIMENT (OSSE) DESIGN 



A. COMPARISONS OF SSM/I IMAGES AND 16 LAYER MODEL OUTPUT 

Before proceeding with the OSSE, EE16 must be shown to 
reasonably simulate SSM/I data. Comparison of model output 
from EE16, the full physics model, with SSM/I observations 
indicates that EE16 can serve as the simulated observation 
system in the OSSE. While EE16 does not perfectly reproduce 
the SSM/I product it does produce storm depictions with two 
useful features. The model output is a good approximation to 
the satellite image when considering storm structure and 
intensity compared to observations. Position and magnitude 
comparisons of post -storm analysis, satellite images and model 
outputs support these findings. 

1. Satellite Data 

Using NRL algorithms and NPS IDEA Lab systems, images 
of SSM/I surface wind speeds and water vapor (total 
precipitable water) are created. SSM/I passes oveir the domain 
are available for 0000 UTC, 0930 UTC and 2200 UTC as shown in 
Figures 4, 5, and 6. With land masking applied to the images, 
the highest concentrations of water vapor, rainfall intensity 
and surface wind maxima are tracked for location, intensity 
and translation speed analysis. These are compared to 
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simulations from the 16 -layer NRL model with a modified 
cumulus parameterization scheme. 

2 . Wind Speed 

SSM/I surface wind speed for 0000 UTC, 0930 UTC and 
2200 UTC 4 JAN 1989 are shown in Figures 8, 9 and 10 
respectively. Unfortunately, most of the wind speed maxima 
indicated in these images are collocated with heavy rain 
cells, reducing the accuracy to within 5 m/s as indicated in 
Figures 4 through 6. In Figures 9 and 10, the wind maxima are 
at the western edge of the swath. In truth, the storm center 
was further west, out of the area of SSM/I coverage. The 0930 
image is straddled temporally by EE16 plots at 0600 UTC and 
1200 UTC, which depict a circular high wind speed area 
propagating northeastward {Figures 11 and 12) . Rainfall near 
the center of this system makes SSM/I magnitudes unreliable, 
but towards the outer edges of the storm there is agreement in 
the contours of 6 to 24 m/s winds. EE16 generated wind fields 
at 1800 UTC and 2400 UTC (Figures 13 and 14) show the same 
comma pattern exhibited in the SSM/I image. (Throughout the 
analysis, "2400 UTC" will refer to the ending time of the 
experiment trials. The actual calendar time would be 0000 UTC 
5 January 1989.) In both model generated plots, the actual 
storm center and associated wind speed maximum lay to the west 
of the comma structure, just outside the satellite swath. 
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Figure 8. SSM/I wind speed (m/s) at 0000 UTC. 



34 





18-24 
12-18 
6-12 
< 6 



Figure 9. SSM/l wind speed (m/s) at 0930 UTC. 
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Figure 10. SSM/I wind speed (m/s) at 2200 UTC. 
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Figure 11. EE16 generated surface wind field (m/s) at 
0600 UTC. 
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Note the wind speed minimum evident at the extreme western 
swath edge near the comma head in the SSM/I image shown in 
Figure 10. This corresponds to a model indicated minimum 
eastward of the comma structure. 

3. Water Vapor (Total Precipitable Water) 

SSM/I water vapor images and total precipitable water 
plots generated by EE16 support the use of EE16 as a simulated 
observational system. The terms 'water vapor' and 'total 
precipitable water' may be used interchangeably in this case. 
Figures 15, 16 and 17 are SSM/I images of water vapor for 
0000, 0930 and 2200 UTC respectively. Figures 18 through 21 
show that EE16 produces a structure very similar to the SSM/I 
images regarding total precipitable water distribution. 
Evolution of the circular precipitable water center defined by 
the 30 -36 kg/m^ contour in Figure 16 to the banded, T-bone 
structure defined by the same contour in Figure 17 compares 
favorably with the EE16 series. In Figure 18, EE16 produces 
a circular 36 kg/m^ precipitable water area centered near 32° 
N, 70° W evolving to a T-bone structure defined by the 30 kg/m^ 
contour in Figure 21. 
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Figure 15. SSM/I water vapor image at 0000 UTC. Contours 
depict precipitable water amounts in kg/m^. 
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Figure 16. As in Figure 15, but at 0930 UTC. 
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Figure 17. As in Figure 15, but at 2200 UTC. 
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at 0600 UTC. 
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Magnitudes of total precipitable water do not support use 
of EE16 as an observation simulator, but the trend in 
development makes them suitable. Since the actual storm 
deepened to 936 mb by 2400 UTC while the 16 layer model goes 
only to 956, it is not expected that all parameters show 
agreement. The SSM/I, with an error of 2.4 kg/m^, indicates 
an initial rise in water vapor followed by a leveling at 45 
kg/m^. The sixteen layer model agrees with this rise in 
precipitable water amounts. Again, these SSM/I values are 
amplified due to the removal of the precipitation screen. 
Since this storm deepened much more than the 16 layer model 
predicts, the SSM/I values are assumed to be much closer to 
the truth. However, EE16 does duplicate the trend and 
approximate the magnitude, and will be accepted as an 
observing system simulator. 

B. EXPERIMENT DESIGN 
1. OSSE 

To observe the impact of satellite data on the NRL 
model forecast, an Observing Systems Simulation Experiment 
(OSSE) is conducted. This OSSE is specifically organized to 
contrast dynamic and diabatic forcing variables. The 
experiment series structure follows the Air Force Geophysics 
Laboratory description of a "Fraternal Twin" OSSE (Lipton, 
1989) . The reference atmosphere will be generated by allowing 
a degraded version of the NRL 10 -layer model (MIO) to run over 
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the 24 hour observ'ation period. Degradation of the model 
results in a greater separation in skill level of the two 
systems being run. In these tests, MIO is hindered by an 
imposed 50% underestimation of horizontal water vapor 
convergence in the Kuo convective parameterization scheme. 
Diffusion is amplified by a factor of four. This inhibits 
latent heat release and, consequently, storm deepening. 
Simulated obseirvations will be generated by the more complex 
16 -layer model (EE16) with full model physics. This version 
models the planetary boundary layer with seven vertical 
levels, compared to MlO's bulk PEL parameterization contained 
in the single lowest level. These simulated observations will 
consist of atmospheric moisture content and surface wind 
speed. 

To conduct the OSSE, these simulated observations are 
then incorporated into the 10 layer model reference atmosphere 
in the various schemes summarized in Table 3 . Each column 
presents a feature of the experiment. Test names start with 
"T" for assimilation of a variable over the entire domain or 
"S" for assimilation only over the areas of satellite coverage 
(swaths) . The second letter identifies the tested parameter, 
as "W" for surface winds, "M" for substitution of the 16 layer 
moisture into the 10 layer model, or "V" for vertical 
distribution of the moisture correction suggested by the 16 
layer model. Column 2 indicates the number of "observation" 
updates made during the experiment. Column 3 states the 
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amount of time each "observation " was held valid. The "Lead" 
column defines the amount of forward lead time allowed in the 
incorporation of 16 layer data. (For example, the 0600 UTC 
EE16 winds are used by MIO at run time 0300 for a lead of 3 
hours . ) 

TABLE 3. SUMMARY OF EXPERIMENTS 



Test 


Observations 


Period 


Area 


Lead 


Ending 

Pressure 


TWl 


4 (Surface 
Wind) 


6 Hours 


Domain 


None 


966 mb 


SWl 


4 (Surface 
Wind) 


6 Hours 


Swath 


None 


969 mb 


TV2 


4 (Vertical 
Moisture) 


6 Hours 


Domain 


6 Hours 


974 mb 


TM2 


4 (Moisture) 


6 Hours 


Domain 


6 Hours 


968 mb 


SM2 


4 (Moisture) 


6 Hours 


Swath 


6 Hours 


973 mb 


TV3 


4 (Vertical 
Moisture) 


6 Hours 


Domain 


3 Hours 


973 mb 


TW4 


2 (Surface 
Wind) 


12 Hours 


Domain 


6 Hours 


967 mb 


TM4 


2 (Moisture) 


12 Hours 


Domain 


6 Hours 


973 mb 


TV4 


2 (Vertical 
Moisture) 


12 Hours 


Domain 


6 Hours 


973 mb 



2. Analysis plan 

Various sensitivity tests will be run using 
combinations of simulated satellite parameters. Although the 
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SSM/I provides only directionless wind speed, since these 
speeds are successfully modeled it will be assumed that the 
model generated wind directions are also valid. Another 
induced artificiality is the vertical distribution of water 
vapor. Since the SSM/I provides a single value for water 
vapor representing the contents of the atmospheric column 
above a point at the surface, a method of numerically 
distributing this water throughout the column is required. 
Section D contains a description of the vertical distribution 
scheme employed by the 10 layer model. 

Following collection of NCAR-Graphics generated plots 
of various fields for each model run, development of the storm 
will be examined to assess the relative importance of wind 
speed and the horizontal and vertical distribution of water 
vapor. Besides the storm features described by other studies 
in the previous section, analysis of the model runs will 
focus on the meteorological processes listed below. 

Vertical storm development will be examined with 
respect to the increase of circulation intensity with height 
and interactions with short waves. Wind structure, pressure 
and temperature contours will be used to determine the 
character of the flow and the structure of the storm. 

Of particular interest in this study is the structure 
and function of moisture distribution within the storm. 
Diurnal effects of moisture transport will be observed; 
however, conclusions will be limited since this is a single 



52 



24 -hour period under discussion. Cross sectional analysis 
will describe vertical motion fields in the vicinity of 
fronts . 

Four series of sensitivity tests are performed on the 
ten layer model (MIO) to determine the most critical 
parameters in describing the behavior of a rapidly deepening 
cyclone. Each series of tests contains up to three variations 
on variables used and spatial domain included. The first 
series updates the MIO fields with coinciding fields from EE16 
every six hours. Series II initializes the MIO run with 
selected fields from the 0600 UTC EE16 run, and continues to 
update every six hours. Series III resembles the first 
series, with the update times moved three hours earlier. 
Finally, Series IV suggests the pattern obtained if updates 
were received from a polar orbiting satellite, using two 
updates spaced twelve hours apart. Figure 22 is a line 
diagram illustrating the flow of the tests. These tests 
provide information on which parameters MIO is most sensitive 
to in explosive cyclone prediction, but they do not represent 
operationally usable methods since 'a priori' knowledge of the 
atmosphere is assumed in most of the runs. 

C. SERIES I; COINCIDENTAL FOUR PERIOD NUDGING 

In the first experiment series, MIO marches through time 
normally with respect to all the state variables except the u 
and V surface wind component fields. These fields are nudged 
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SENSITIVITY TEST SUMMARY 
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Figure 22. Sensitivity test summary. Asterisks indicate update 
points using EE16 data. Numbers below asterisks indicate 
generation time of update information by EE16. 



towards their respective EE16 -generated values at 0600 UTC 
during the first six hours of the model run, toward the 1200 
UTC value during the second six hours, the 1800 UTC fields 
during the third four hour period, and towards the 2400 UTC 
values for the final four hours of the test. This series 
includes two runs, TWl and SWl. In TWl, winds are nudged over 
the entire spatial domain of the case study, while SWl nudges 
the winds only over the SSM/I swaths covering the ERICA area. 

The primary step in the nudging process involves 
converting the 16 layer wind into a compatible 10 -layer 
product from which a surface value can be extracted. The 
subroutine MWWIND in the NRL Model applies a mass -weighting 
function to the 16 layer vertical wind profile. Figure 23 
illustrates this vertical weighting function. The lowest 
seven layers of the weighted profile are averaged and 
extracted as the "Observed" surface wind. Once this value has 
been obtained, standard Newtonian nudging practices are 
followed to gradually update the MIO winds. In practice, the 
new MIO wind depends on the horizontal pressure gradient, 
Coriolis parameter and "observed" wind. 

For SWl, the model defines the slope of the swath 
boundaries. At the same times as in run TWl, winds over the 
swath region are nudged to coincident EE16 values. Estimation 
to the actual swaths are made, using the 0000 UTC swath for 
model nudging from 0000 to 0600, the 0930 swath for data 
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16 Layer data to 10 Layers. 
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between 0600 and 1200 and the 2200 swath for requirements 
after 1800. 

D. SERIES II: MOISTURE ASSIGNMENT AND DISTRIBUTION 

The second series of experiments diagrammed in Figure 22 
begins with assignment of the 0600 EE16 moisture ("q" values) 
to the 0000 UTC MIO atmosphere. The model runs until 0600 
when the 1200 EE16 moisture is similarly assigned. This six 
hour lead repeats until the end of the experiment. Both a 
total domain test (TM2) and swath assignment test (SM2) are 
performed. 

Series II also includes a procedure for distributing the 
total precipitable water "observed” over the vertical 10 layer 
column (TV2) . The NRL Model subroutine QDISTNRL performs the 
vertical moisture distribution adjustment at those grid points 
indicating a significant difference between MIO and EE16 
derived total integrated water content. This correction can 
be negative to indicate drying of the column, or positive for 
moistening. At these points, vertical levels in the 10 layer 
model are examined to determine unsaturated layers which could 
absorb portions of a positive (adding moisture) correction, 
using subroutine ESAT to provide the specific humidity profile 
at the grid point. The MIO value for humidity in the column 
is then adjusted by the following: 
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Where q(K)^^^ represents the new moisture value in a given 
layer, CORR the total correction applied based on the EE16-M10 
moisture difference, RELH the relative humidity, and RELH^un^^t 
the relative humidity in the unsaturated layers. After the 
vertical moisture redistribution, control checks insure no 
layer has been supersaturated or "hyper-dried”. If such has 
occurred, humidity profiles are recomputed using the 
adjustment scheme until no supersaturated or "hyper-dried" 
regions exist. 

E. REMAINING SERIES 

1. Centered Four- Period Assignments 

Series III consists of a single vertical distribution 
run, TV3 , using the algorithm described in QDISTNRL. Here, MIO 
runs normally until 0300 when the integrated water content 
computed by EE16 at 0600 is inserted. MIO runs from this 
assignment until 0900, reacting to the redistribution. Using 
this scheme, the observation is centered within the 
assimilation period (i.e., the 0600 observation is centered 
between 0300 and 0900) . Three periods of six hour run times 
follow, centering the 1200, 1800 and 2400 EE16 integrated 
water content observations as shown in Figure 22. 
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2. "Polar Orbiter" Simulator 

Series IV contains three runs over the entire domain, 
testing winds (TW4) , surface moisture (TM4) and vertical 
distribution moisture adjustment (TV4) when "updated" twice at 
a twelve hour interval as opposed to the six hour updates of 
the previous tests. A centering scheme is used again, as MIO 
is assigned the 1200 UTC EE16 value at model time 0600 UTC, 
and assigned the 2400 UTC EE16 value at 1800 UTC, allowing the 
model to run normally for the final six hours. 

F. FINAL COMMENTS 

The experiments described in this chapter are not 
operationally feasible since most require assimilation of 
"future" data which would not be available to an operational 
model. They can, however, provide insight toward the trend of 
corrections needed to improve the model performance for a 
given type of storm. In the next chapter, these sensitivity 
test will show the vulnerable areas of MIO in modeling 
explosive extratropical cyclones. 
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V. ANALYSIS OF RESULTS 



A. OVERVIEW 

ERICA IOP-4 exemplifies a dynamically -driven storm with 
diabatic processes playing a secondary role in the development 
of the storm. The assimilation of wind speed data has the 
greatest effect on model performance. Model depiction of 
cyclone circulation interacting with the environmental flow 
impacts the development rate of the storm. Secondary 
deepening effects are revealed in precipitation amounts and 
areal coverage which imply the degree of latent heat release. 
Although no experiment deepens the storm to the analyzed 936 
mb intensity, or even to the 956 mb low predicted by the EE16 
run, significant improvements from the 971 mb degraded run are 
observed with wind speed inputs. Greater deepening rates with 
better defined frontal structures are gained from wind speed 
inputs. Sensitivity to moisture distribution is also 
demonstrated. 

B. COMPARISON OF SENSITIVITY TEST RESULTS 

Four areas of study arise from the sensitivity tests 
described in Figure 22. First/ although modelers acknowledge 
the value of accurately initialized winds in a numerical 
prediction scheme, these values are not often readily 
available. Sensitivity test Series I (TWl vs. SWl) examines 
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the degree of improvement afforded by incorporating the winds 
over the satellite swath compared to the improvement gained by 
incorporating winds over the entire domain. Second, 

sensitivity test Series II provides two comparisons to 
determine the most productive use of satellite obtained 
moisture. Within this series, a comparison of swath vs. 
domain surface "q" values is performed similar to the first 
wind test. Additionally, the series tests the value of 
vertically distributing the "observed" moisture throughout the 
atmospheric column in a domain test. 

Test Series III, consisting of one vertically distributed 
moisture test, compares to the vertical test from Series II 
giving an indication of the model sensitivity to assimilation 
time difference. This comparison investigates the value of 
reducing the six hour lead time to three hours as well as a 
"centered" estimation time versus a "lead" time. For the 
final comparison. Series IV places surface winds, surface 
moisture and vertically distributed moisture in identical 
update situations to determine the most effective of the 
three parameters for explosive cyclone modeling. 

C. SERIES I: SWATH VS. DOMAIN WIND ASSIMILATION 

Storm development is traced for the first series of 
experiments in Figure 24. Significant differences in wind 
speed structure develop at 1200 UTC. At the surface in the 
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Figure 24. Series I storm development. 
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immediate vicinity of the storm (directly east of the DELMARVA 
peninsula at 38° N, 64° W) , TWl models a strong zonal wind 
speed gradient, shown in Figure 25. At the same time, SWl 
displays more curvature as the local wind speed high of 32.2 
m/s has already swung to the southeast of the wind speed 
minimum. Winds above 500 mb show little change and do not 
appear to be sensitive to the surface wind distribution this 
early in the test. 

The experiments in Series I provided the best storm 
development in terms of surface pressure deepening. Both runs 
have a significant impact upon the surface momentum field. 
Figure 26 presents surface and temperature fields at 2400 UTC. 
With swath wind inputs from 2400 UTC applied, SWl only lags 
TWl by 3 mb in storm development. SWl's ability to remain 
close to TWl in storm development is largely attributable to 
the early swaths coinciding with the areas of most intense 
storm development. Although TWl surface pressure remained 
consistently deeper than SWl, the latter displays more 
distinct frontal structure. Since the first satellite pass 
covers the area of storm development, the cyclone is 
highlighted by the model in SWl. The two experiments produce 
wind fields with differing orientations. Figure 27 displays 
the surface isotach field for each run at 2400 UTC. TWl 
places local maxima of 34.5 m/s and 35.0 m/s north and south, 
respectively, of the local minimum. SWl positions local 
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Figure 26. TWl/SWl surface pressure (solid lines, mb) and 
temperature (dashed, °C) at 2400 UTC. 
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Figure 27. TWl/SWl isotachs (m/s) at 2400 UTC. 
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maxima of 35.1 m/s (east) and 33.7 m/s (west) around the local 
minimum. Comparison of precipitation areal coverage and 
amount shows small differences related to LHR and upward 
motion between the two experiments. Although a specific study 
of heating rate assimilations is to be done in a future paper, 
impact of differing heating rates is evident in these 
experiments . The improved deepening rate of TWl can be traced 
back to the 1200 UTC cross sections through the warm front. 
In Figure 28, on the southern edge of the warm front the omega 
field shows upward motion and tighter omega gradients to 
higher levels than the SWl plots shown in Figure 29. The 
tighter gradient feature continues through 2400 UTC, ending up 
farther north in TWl than SWl. Across the cold front, which 
is better represented in SWl due to the contrast in resolution 
inside and outside the swath, rising motion in the warm sector 
at 12 00 UTC is distinct in SWl and minimal in TWl. This 
contrast does not maintain itself past 1500 UTC, as the 
subsequent satellite pass misses portions of the storm as the 
low center moves eastward. 

Statistical analysis is performed to check the performance 
of the model and the impact of the experiment variations on a 
grid scale over the entire domain. Statistical plots used in 
this analysis consist of root -mean- square (RMS) errors and SI 
scores computed against the corresponding EE16 fields. In 
Series I, Figure 30 shows surface wind RMS errors are 
initially larger for TWl than for SWl. By the end of the run. 
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Figure 28. Warm front cross section of omega 
(microbars/sec) for TWl at 2400. "A" and "B" correspond to 
points on Figure 26. 
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it is evident that TWl has maintained an overall lower error 
compared to the EE16 run. During the course of the 
experiment, the swath passes used in SWl get progressively 
farther from the storm center and the highest winds. By the 
end of the experiment, the swath update information is 
correcting low magnitude wind differences while leaving the 
major discrepancies untreated. At 850 mb, both runs have 
similar RMS wind errors indicating the relatively low impact 
the domain- swath variation has on vertical wind structure. 
Finally, Figure 31 shows the comparison of storm tracks based 
on low pressure centers generated by the two experiments. SWl 
follows the same track as TWl but at a slightly reduced speed. 
Again, this is swath location dependent as the earlier updates 
improve the storm center information, while later updates 
outside the cyclone reduce the environmental wind speed 
gradient . 

D. MOISTURE SENSITIVITY 

1. Domain vs. Swath Comparison 

From Table 3 in Chapter IV it is clear that 
incorporation of surface moisture over the entire domain is an 
improvement from incorporation of swath moisture alone. 
Experiment TM2 produces a final pressure of 968 mb compared to 
SM2's 973 mb. Moisture difference plots show that both runs 
produce changes of equal magnitudes. The structure, however, 
differs significantly. Figure 32 compares the precipitable 
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RMS Surface Wind: Series I 
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Figure 30. Series I surface wind RMS errors. 
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Figure 31. SWl vs. TWl storm tracks. 



water fields for TM2 and SM2 . TM2 displays a broken 
structure behind the cold front in the area of expected cold 
air advection. This may be an attempt by the model to display 
sub- grid scale convection. SM2 maintains a smoother 
integrated water pattern in the same region. The central 
surface pressure of the storm throughout these two experiments 
remains close, within 2 mb as shown in Figure 33, until the 
final three hours when TM2 drops 7 mb compared to SM2"s 2 mb. 
Differences in storm structure leading up to this disparity 
become evident at 1800 UTC when 2400 UTC moisture values are 
assigned. In the swath test. Figure 34 shows temperature 
contours near the cold front (trailing south of the main low 
following the 18° and 20° C isotherms) are more closely packed 
in the SM2 test. The cold advection area behind the front is 
very jagged and disorganized. At 2100 UTC, the irregular 
temperature pattern in this cold advection remains, and 
development of the storm appears stalled. Figure 35 displays 
the large areal precipitation which may be associated with 
small scale convection in TM2 . 

Wind speed statistics show a marked difference between 
the two runs in the 500 mb wind speed field for the 1800 to 
2400 UTC period. At 1800 UTC, when the 2400 UTC moisture is 
assigned. Figure 36 shows the a jump in RMS wind speed error 
for SM2 that ends up exceeding the error in the control run. 
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Figure 32 . Precipitable water content (xlO kg/m^) for TM2 
and SM2 at 2400 UTC. 
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Figure 33. Series II storm development. 
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Figure 34. SM2/TM2 surface pressure (solid, mb) 
temperature (dashed, °C) at 2400 UTC. 
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and 




TM2 . Contour interval is 10 cm. 



77 



RMS 500 mb Winds: Series II 




- TV2 TM2 « SM2 



Figure 36. Series II 500 mb wind speed error. 
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Increasing RMS errors in 850 mb temperature also begin at 1800 
UTC, but these errors do not exceed those of the control run. 
Finally, Figure 37 shows the 500 mb moisture fields are nearly 
a perfect match throughout the run in TM2. SM2, however, has 
errors nearly an order of magnitude greater, peaking at 1800 
UTC and decreasing. 

2. Surface Only vs. Vertical Moisture Distribution 

Employing the vertical moisture distribution scheme 
described in Chapter IV slows development of the storm 
providing a poorer forecast. Compared to the control run, use 
of the vertical distribution scheme reduced the effectiveness 
of the moisture inputs. Both frontal moisture band structure 
and cold air convection are suppressed. The vertical 
distribution technique improves the temperature field at 850 
mb, but not far above or below. Differences in the surface 
pressure start at 1 mb at 0900 UTC (Figure 33) . During the 
final three hours of the experiment, the surface only test 
(TM2) drops 9 mb, while inclusion of the vertically 
distributed moisture (TV2) nets only a 3 mb drop. TV2 
displays the same unsettled isotherm pattern in the cold 
advection region seen in earlier tests. In integrated water 
vapor comparisons, TV2 shows some of the mottled moisture 
cells, certainly a more accurate depiction than the swath test 
but not as accurate as TM2 . TM2 suppresses the rainfall. 
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RMS 500 mb Moisture: Series II 
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Figure 37. Series II 500 mb moisture RMS error. 
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reducing latent heat release and slowing development. 
Vertical distribution of the moisture tends to slow winds in 
the lower portion of the atmosphere, notably the 34 vs 30 m/s 
winds maxima at 900 mb, 1800 UTC in Figure 38. The RMS error 
for 850 mb winds of TV2 was the worst among the 9 runs, while 
TM2 produced the smallest errors. 

Upper atmosphere temperature is sensitive to the two 
differing methods of determining moisture throughout the 
column. TV2 maintains RMS errors of temperature at 500 mb 
greater than the control run error (Figure not shown.) TM2 
errors match the control run until 1800 UTC, and then improve 
dramatically. Similarly, mid- level moisture errors are 
magnified with the vertical distribution. At 850 mb the RMS 
"q" error for TV2 is worse than the control run until 1500 
UTC. 

As discussed in the description of the vertical 
distribution scheme in Chapter IV, the vertical moisture 
distribution scheme computes the difference in total 
integrated water determined by the 10 layer model and by the 
16 layer model. After the individual layers in MIO are 
checked for saturation, the computed difference is spread over 
the unsaturated layers, each layer receiving a piece of the 
correction proportional to the layer's relative humidity for 
a moistening correction. Conversely, for a drying correction 
the amount of moisture removed from each layer is proportional 
to its relative humidity. The model is checked to prevent 
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removal of moisture in a layer resulting in a negative value. 
This scheme has the effect of restricting the moisture changes 
to the lowest layers (700 mb and below) of the atmosphere 
where the relative humidity is highest. Confining the mass 
changes to these lowest levels seems to have suppressed storm 
development by inhibiting mass and momentum transport. 

E. TIME SENSITIVITY 

No significant surface pressure decrease differences exist 
between TV2 and TVS. Structurally, a better defined cold 
front exists in TVS at 1200 UTC, but the time shift does not 
reduce the structural changes in the cold advection region. 
No significant changes in moisture or rainfall coverage exist. 
In most of the statistical checks, TVS matches the trends of 
TV2 with slightly improved errors. Notable exceptions to this 
include the RMS wind error at the surface. The reduced lead 
time experiment, TVS, again matches the trend of the long lead 
TV2 run, but the error is greater for the shorter lead time. 

With the reduced lead time, vorticity above the boundary 
layer associated with the surface cyclone endures longer. 
Vorticity at 700 mb in TVS goes from 79.SxlO‘^s‘* at 1800 UTC to 
72.2xl0'^s'^ at 2400 (Figure S9), while the drop is larger, from 
77.7xl0'^s'^ to 66.9xl0'^s*^ in TV2 (Figure 40). The decreasing 
of above -boundary layer support contributes to the poor 
performance of TV2 . 
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Figure 40. As Figure 39, but for TV2 . 
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F. INDIVIDUAL PARAMETER COMPARISONS IN THE POLAR ORBITER 



SIMULATION 

Experiment Series IV compares total domain applications of 
each tested parameter. The sensitivity of the model to 
surface wind speeds is readily apparent in this series. TW4 
produces a 967 mb storm, the second best performance of the 
entire experiment group. Experiments TM4 and TV4 each produce 
a 973 mb storm, the second weakest forecast of sea level 
pressure in the field. Figure 41 shows the development of 
these three runs. TW4 improves upon the other runs almost 
immediately, dropping to 3 mb below the moisture experiments 
as early as 0600 UTC. At 1200 UTC development of TW4 slows, 
allowing TM4 to surpass the wind test by 1 mb. Here, 
development of TM4 quickly slows, while the other experiments 
continue to deepen. 

The surface pressure and temperature plot of TW4 in Figure 
42 at 0600 shows the increased area of cold advection around 
the western edge of the low. Good advection can be inferred 
from the north end of the low westward around to the southern 
edge, while in the other two experiments advection extends 
only to the south- southwestern point. At 1500 UTC (Figure not 
shown) , determination of the cold front becomes difficult in 
TW4, while the 16-18 ° C isotherm packing in TM4 and TV4 makes 
the front obvious. Warm advection across the front is 
increased in the moisture runs as well. This suggests, since 
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Figure 41. Series IV storm development. 
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the last update by "observed" fields was with the 1200 UTC 
fields inserted at 0600, that the model reacts quickly to wind 
speed distribution and slowly (here about 9 hours) to moisture 
inputs. Warm advection is strongest in the vertically 
distributed run, TV4 . 

By 1800 warm air advection modeling by the TW4 run 
surpasses that of the moisture tests, creating a marked "S" 
shape baroclinic zone. Significant warm air advection occurs 
northeast of the low, and strong cold air advection occurs 
west -southwest of the low (Figure not shown) . By the end of 
the run, TW4 has produced a storm farther north with stronger 
warm and cold advection, 6 mb deeper than the moisture runs. 

Water vapor distribution remains similar among the three 
experiments until 0900 UTC. Both moisture tests show the 
previous patchy moisture pattern in the cold advection region. 
The wind test remains smooth for the first half of the 
experiment. The contours become disturbed but never fully 
achieve the mottled pattern. Total precipitation in the 
surface moisture run exceeds that of the vertically 
distributed moisture on the surface wind test. 

Incorporation of the 1200 UTC observed surface winds at 
0600 distorts the isotherm field in TW4 (Figure not shown) . 
At this early stage of the storm, two nearly equal wind maxima 
are placed northwest and southeast of the expected storm 
position. The moisture experiments each place a 31.7 m/s wind 
max off Cape Hatteras. This distortion does not extend far 
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into the upper atmosphere. By 1200 UTC, the surface pattern 
among the three experiments are in agreement, although the 
vertically distributed moisture experiment produces slightly 
decreased wind speeds. From this point, TW4 proceeds to 
develop intense speed gradients around the storm, while TV4 
continues to lag TM4 in magnitude. 

Figures 43 through 45 show the 2400 UTC surface pressure 
and temperature fields in these three experiments. Cross 
sections (not shown) taken perpendicular to the warm and cold 
fronts reveal the vertical motion gradients present in the 
baroclinic zones. TW4 has a single intense upward motion 
area, while the moisture runs show multiple weaker upward 
motion cells. Along the warm front, however, TM4 displays the 
most intense vertical motion. 

Statistics of the 850 mb temperature field indicate that 
at this level the moisture runs significantly improve model 
performance (Figure 46) . 850 mb RMS temperature errors in the 
wind test are identical to the control run errors. At 500 mb, 
the vertically distributed moisture run quickly degrades, 
ending up with the largest temperature errors of the field. 
Again, the wind run maintains errors close to the control run. 
Trends of moisture errors are very erratic. While the wind 
and vertically distributed moisture parallel the sawtooth 
control pattern in Figure 47, TV4 degrades quickly to the 1800 
UTC mark, improving slightly by 2400 UTC. At 850 mb, TM4 
stands out as the errant run, this time showing improvement 
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Figure 45. TV4 surface pressure (solid) and temperature 
(dashed), 2400 UTC. 
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850mb RMS Temperature Errors: Series IV 
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Figure 46. 



Series IV 850 mb temperature RMS errors. 
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500 mb Moisture errors: Series IV 




-m- TM4 -■+- TV4 « TW4 



Figure 47. Series IV 500 mb moisture RMS errors. 
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over the control run but without a determinable trend (Figure 
48) . In the 500 mb level, both moisture runs spike at 1800 
UTC, with TM4 improving dramatically afterwards. Finally, RMS 
wind errors (Figure 49) again show a 1800 UTC error peak for 
TM4 at the surface, with improvement by TW4 at 500 mb. 
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RMS 850mb Moisture: Series IV 
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Figure 48. Series IV 850 mb moisture RMS errors. 
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RMS 500 mb Winds: Series IV 
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Figure 49. Series IV 500 mb wind RMS errors. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

Without question, accurate surface wind estimates enable 
the NRL model to predict more accurate, deeper storms. 
Generally, total domain coverage with observed winds is 
desired, but swath winds provide significant improvements to 
the forecast if the swath falls over the area of main storm 
activity. From this information, frontal structure and 
temperature advection in the most important energy transfer 
regions can be modeled. 

Incorporation of observed moisture as done in this study 
fails to consistently improve the prediction. Using relative 
humidity as a distribution variable places most of the 
moisture correction at the lowest levels of the atmosphere. 
Placement of too much moisture at the wrong levels inhibits 
momentum transfer and inhibits the storm development. Upper 
level winds decrease, reducing horizontal advection and 
vertical motion. The advantage of observed moisture is the 
improvement in open cell convection modeling as seen in the 
water vapor distribution fields. 

Specifically, the model exhibits sensitivity to all tested 
parameters. Use of moisture updates enhance the structural 
modeling ability of MIO, especially behind the cold front. 
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The magnitude of these corrections compared to the control run 
are small, not approaching the values given by the satellite 
images as quickly as hoped. Still, MIO itself does not 
produce cellular moisture structures behind the cold front 
without satellite information. Vertically distributing the 
moisture reduces the accuracy of the surface water vapor 
field, but improves the midlevel temperature depiction by 
manually redistributing the heat source poorly modeled by MIO. 

Assimilation of winds improves the momentum transfer from 
the lower levels, which has the greatest effect on the 
development of the storm. This carries to the upper 
atmosphere, as shown by the substantial improvement in 500 mb 
vorticity. Sensitivity to the surface winds is enhanced by 
strong vertical coupling of the momentum fields in MIO. This 
coupling does not occur instantly. Based on the earlier times 
of the swath- domain wind experiments in Series I, the model 
requires about six hours to assimilate information over the 
domain, while local effects are evident in as little as three 
hours . 

MIO is highly sensitive to the placement of observed data. 
Wind and moisture information from areas of central storm 
development improve the forecast. Swaths which cover more 
environmentally flow- controlled areas inhibit storm 
development by reducing gradients from storm center to the 
outlying environment. 
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The NRL model exhibits strong sensitivity to wind speeds. 
In this dynamic case, updating the surface momentum field 
carries over to heat transfer and mass continuity processes. 
Weaker sensitivity to moisture, which also affects heat and 
mass distribution, suggests that momentum fields are the most 
useful parameter in the modeling of the explosive cyclogenesis 
in this case. 

B. RECOMMENDATIONS 

Future studies may enjoy the opportunity of directly 
implanting satellite derived data into the reference 
atmosphere model instead of performing an OSSE. In that 
scenario, wind speeds obtained using current algorithms may be 
weighted according to the wind speed accuracies presented in 
Chapter II. This will allow the model to lend more credence 
to the most accurate winds. An algorithm or parameterization 
scheme may be devised to incorporate wind directions as well. 

Certainly, improvements to the vertical distribution 
scheme are needed. Relative humidity seems to be a poor 
indicator to use in deciding where to place additional 
moisture. Other candidates include saturation mixing ratio or 
potential temperature. When later versions of the DMSP 
satellite become operational, an on-board microwave profiler 
will provide actual vertical moisture profiles. 
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